Acoustic attenuation rate in the Fermi-Bose model with a finite-range 

fermion-fermion interaction 
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We study the acoustic attenuation rate in the Fermi-Bose model describing a mixtures of bosonic 
and fermionic atom gases. We demonstrate the dramatic change of the acoustic attenuation rate as 
the fermionic component is evolved through the BEC-BCS crossover, in the context of a mean-field 
model applied to a finite-range fermion-fermion interaction at zero temperature, such as discussed 
previously by M.M. Parish et al. [Phys. Rev. B 71, 064513 (2005)] and B. Mihaila et al. [Phys. 
Rev. Lett. 95, 090402 (2005)]. The shape of the acoustic attenuation rate as a function of the boson 
energy represents a signature for superfluidity in the fermionic component. 
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I. INTRODUCTION 



Recent realization of superfluidity in a dilute gas of 
ultracold fermionic atoms, is driving renewed theoretical 
efforts aiming at a quantitative understanding of many- 
body correlations in a dilute Fermi gas. At low density, 
for finite-range interactions, one expects that all sensitiv- 
ity to the detailed features of the interaction is lost 
and the system's quantum behavior depends only on the 
scattering length and the range of the interaction. Fur- 
thermore, as we approach the "unitarity" limit @, i.e. 
the limit where the scattering length is much larger than 
the mean inter-particle separation, the system exhibits 
universal behavior and the energy is determined entirely 
by the density, leading to high sensitivity to many-body 
correlation effects. 

Built on the successful realization of quantum degen- 
erate systems consisting of boson and/or fermion Q 
atoms, experimental progress in cold atom physics has 
recently resulted in the generation of ultracold molecules 
from ultracold atoms [f|, Bose- Einstein condensates 
(BEC) of molecules @, and the first experimental re- 
alization of the small pair to large pair crossover @, H, Q • 
In one type of experiments @, S, @] , the system consists 
of one-type of cold, degenerate, fermionic atoms in an 
optical trap. The strength of the interatomic interaction 
can be continuously adjusted via a magnetically-tuned 



Fcshbach resonance between a closed and an open chan- 
nel. Hence, the system can evolve from a BEC phase in 
which spatially non-overlapping (Shafroth) pairs [l(| are 
bound together, to a Bardeen-Coopcr-Schricffcr (BCS)- 
type superfluid [ll[ involving correlated atom pairs in 
momentum space. 

Significant experimental work has been also directed 
towards the study of mixtures of bosonic and fermionic 
atom gases [H, [l3| . Such mixtures remind of the strongly 
correlated 3 He- 4 He mixtures [l4[, and can be obtained by 
sympathetically cooling a gas of fermionic atoms by using 
a condensed boson g ClrS CIS 0/ refrigerator. In this paper, we 
will focuss on the latter type of system, and discuss the 
fermionic response to a density fluctuation of the Bose 
condensate in the trap, with special emphasis on the ef- 
fects due to the finite-range character of the fermion- 
fermion interaction. We will argue that the shape of the 
acoustic attenuation rate depends on the BEC or BCS 
character of the fermionic liquid. 



II. ACOUSTIC ATTENUATION IN THE 
FERMI-BOSE MODEL 

The Hamiltonian describing the interacting Fermi- 
Bose system is given as the sum of a boson, fermion and 
interaction components 
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This model, proposed in this context by Timmermans et 
al (l5j . is similar to the a model introduced by Bardeen 
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et al [l6| to describe the effect of the fluctuations of the 
polarization in an excitonic liquid and a gas of electrons. 
This model has since been the subject of intense investi- 
gations in connection with the Feshbach resonance mech- 
anism in ultracold fermionic gases [l?], G3 ■ 

In the weak coupling limit, the Fermi and Bose compo- 
nents can be treated independently, and the many-body 
wave function of the Fermi-Bose mixture is the direct 
product of independent Fermi and Bose wave functions, 
\®fb) = \®b, 3>f)- The dispersion relation of the dilute 
Bose condensate is [ljl 



(2) 



where c = \J \b Pb /ms is the phonon (BEC excitation) 
velocity of sound, = 1/ 'y 1 AXePb^b is the boson co- 
herence length, and Xb = Anas/mB is the strength of 
the bosonic pseudo-potential. The coherence length, i.e. 
the length scale at which the condensate varies spatially, 
is typically of the order of £b ~ 0.1 — 1 /mi, while c is 
approximately 0.01 /imMHz in condensates with density 
of 100 fim~ 3 . The scattering length, as, and the den- 
sity of the bosonic component, ps, arc tunable param- 
eters which can be used to modify the dispersion rela- 
tion ©, but c£ B = 1/(2 ma) j c/£ B = 2 Pb ~ Pb a B 
and, therefore, the parameters c and £g cannot be varied 
independently. Throughout this paper H = 1, i.e. the 
energy is measured in units of frequency. 

The low-energy state of the system is described by the 
free-energy associated with the BCS-reduced Hamilto- 



H F = ^2 £ k (a kT akt + a k;i a ki) 
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with ek = k 2 /(2mF). Then, in a zero-temperature mean- 
field approximation [2(| the fermionic ground- 
state in the Hartree-Fock-Bogoliubov representation [23[ 
is the standard BCS variational wave function given by 



|<M = TV exp S~] — 




(4) 



The BCS ansatz interpolates smoothly between the BCS 
and BEC limit. The parameters {«k, v^} are obtained 
by solving the system of equations = \ — (ek — 
(j,)/(2Ek) and Ukfk = A^/(2E-\ S ), subject to the nor- 



malization condition, lukl 



«k 



1. Here, /i is the 



fermionic chemical potential such that spin-up and spin- 
down states are equally occupied, and denotes the 
quasi-particle spectrum in the fermionic ground state, 
= (ek — m) 2 + A k - The pairing gap, Ak, is obtained 
as A k = -EpVi-pAp/fip. 

The interaction between the Fermi and Bose compo- 
nents is assumed to be weak, and the interaction Hamil- 
tonian, H; nt , describing the interaction of a single acous- 
tic phonon of momentum q and energy uj q with the 



fermionic component can be written similarly to the 
electron-phonon Hamiltonian [24|, as: 



H; 



2p 



(5) 

where p is the mass density of the composite gas. (For 
an alternate approach based on a density-density inter- 
action, see Ref. [25|.) In first-order perturbation the- 
ory |26[ , the acoustic attenuation rate is given by 



T(q) = 2tt |(i|H int |/)| 2 S(Ei-E f ) 



(6) 



At zero temperature, for a spin-unpolarized Fermi gas, 
we obtain 12411: 



IjJq 



(7) 



q 2 f d 3 k 

with k' = q + k. In the mean-field approximation, the 
normal and anomalous densities in the fermionic ground 
state are defined as pk = (^| a kt ak Tl^) = l u k| 2 , and ^k = 
($|a_kj.ak| |$) = w k u k- The acoustic attenuation rate 
appears as the superposition of components 7 P and 
which depend separately on the normal and anomalous 
densities, respectively (see also Ref. p7|). 

In principle, one should also consider the possibility 
of exciting the Andcrson-Goldstonc mode of the fcrmion 
superfluid system. However, as it has been discussed in 
Ref. (25[, the Anderson- mode excitation can be avoided 
by specifying that the BEC velocity of sound, c, is less 
than the velocity of the Anderson mode, v, i.e. we require 
c < v. Then, the excitation of the collective Anderson 
mode by a BEC phonon mode is energetically forbid- 
den. This is a reasonable assumption for the mixtures 
of bosonic BEC and fermionic atom gases under con- 
sideration here: In the BCS regime, the velocity of the 
Anderson mode, v, is equal to the Fermi- velocity vf and, 
for typical atom trap conditions, that can be easily sat- 
isfied. The only cold atom mixtures for which c > Vp 
contain ultra-low density fermion systems or ultra-high 
density BEC's: the first type would be difficult to re- 
alize because of the very low critical temperature for 
fermion superfluidity, the second type would be difficult 
to achieve experimentally because three-body recombi- 
nation processes would quickly deplete the experimental 
system. Hence, while the Anderson mode will most likely 
be somewhat softened in the crossover regime (since we 
expect it to go over into the Bogoliubov mode of the 
BEC-system in the BEC limit of the crossover), by choos- 
ing a sufficiently weakly interacting BEC or low density 
BEC, it should always be possible to ensure that c < v 
throughout the entire crossover. 
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FIG. 1: (Color online) Acoustic attenuation rate components 
corresponding to the toy model described by Eq. (JSJ . Here we 
study a large pairing-gap, Ao/sf, regime: we set Ao/sf = 
4, fi/eF = 1, and ^b^f = 5. 



FIG. 2: (Color online) Acoustic attenuation rate components 
corresponding to the toy model described by Eq. (JSj . Here 
we study a small pairing-gap, Aq/ef, regime (BCS-like): we 
set Aq/ef = 0.1, h/ef = 1, and £sfcF = 5. 



III. A SIMPLE MODEL 

In order to study the phase-space constraints subject- 
ing the acoustic attenuation rate, we consider first a 
schematic model of the momentum-dependence of the 
pairing gap, suggested by Comte and Nozieres [2ll |: 



A* 



A 



i + (k/k F y 



(8) 



Note that the realistic momentum dependence of the 
pairing gap, as predicted by the one-channel fermion- 
fermion model with a finite-range interaction to be dis- 
cussed in the next section, has a much lon ger tail than 
the simple ansatz of Comte and Nozieres [21|. Unless 
otherwise specified, in the following we set Hckp/sF = 1 
for simplicity. We also keep fixed the chemical potential, 
h/ef, and the coherence length of the Bose-gas compo- 
nent, £,BkF, and vary the pairing-gap parameter Aq/ef- 
For illustrative purposes, in Figs. [T] and [2J we plot the 
energy-dependence of the various attenuation-rate com- 
ponents, 7 K and j p , for the regime of large and small 
values of the pairing-gap, Ao/e_f, respectively. Because 
the quasi-particle spectrum in the model described by 



3 




Eq. 



has a minimum at finite momentum, k m , then 



FIG. 3: (Color online) Shapes of the energy-dependence of the 
attenuation-rate as a function of the pairing-gap parameter, 
Ao/ef, for fixed values of /j,/ef = 1 and £sfcF = 5. 



the conservation of energy, Ey = uj q — E^, and momen- 
tum, k' = q + k, leads to four branches. The available 
phase-space satisfying these constraints leads to kinks ob- 
served in the energy-dependence of the acoustic attenu- 
ation rate. In particular, we have four "critical" points, 
as depicted by the vertical lines in Figs. [T] and [21 In 
the increasing order of the "critical" energy, we have: 
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u) q< i the energy value at which the momentum conserva- 
tion corresponding to the second (k > k m , k' < k m ) and 
third (k < k m , k' > k m ) branches, is no longer satisfied. 
This failure generates the most remarkable feature in the 
acoustic attenuation rate, resulting in a decrease in the 
acoustic attenuation rate; ujq^ the energy value at which 
the fourth branch (k < k m , k' < k m ) also fails due to 
failed momentum conservation; uj qt s the energy value at 
which the fourth branch vanishes; and 0J q ^ the energy 
value at which the second and third branches also van- 
ish: from this point on, the only contribution is from the 
first branch (k > k m , k' > k m ). 

The energy distribution of the acoustic attenuation 
rate may feature a "shoulder" just before 0^,3, because 
the domain of the integrals corresponding to branches 
three and four is k G [0, k m ], i.e. the integration domain 
for these branches is cj g -indcpendent (see circled area in 
the top panel of Fig. [TJ This "shoulder" is apparent in 
the large Aq/ep regime, and disappears as we approach 
the normal Fermi gas limit. 

In the small Aq/ep regime, the "critical" energies uj q .\ 
and Ug t 2 arc close together, while the "critical" ener- 
gies io q ^ and w 9i 4 are pushed to relatively higher values. 
When Ao /ef becomes very small, only one of the eight 
contributions to V is significant (due to the fact that the 
anomalous density is a delta-like distribution centered 
around kp), and the kinks in the energy distribution of 
the acoustic attenuation rate are entirely obscured, ex- 
cept for UJq.l- 

In Fig. [3] we show evolution of the attenuation rate 
profile as a function of the pairing-gap, Ao/sp, at fixed 
chemical potential. The shape of the attenuation rate 
changes smoothly between low and large values of the 
pairing-gap. The signature of the crossover is the change 
in the sign of the slope of the acoustic attenuation 
rate, for phonon energies between the threshold energy, 
u) q ^T = 2.Efc m ) and the first "critical" value, uj q y. the 
slope is positive in the small Aq /sp regime, and becomes 
negative in the large Aq/ep regime. From Figs. [1] and [2] 
follows that , the crossover is characterized by the ratio 
of 7k /7p, i.e. the ratio of the anomalous- to the normal- 
density. In the small Aq/ep limit this ratio is small (it 
vanishes for the Fermi gas). In the large Ao/sp limit, 
the anomalous-density contribution is comparable in size 
(but still smaller) than the normal-density contribution. 
In the limit when Ao — > 0, we recover the normal Fermi 
gas result. 



IV. EFFECTS OF A FINITE-RANGE 
FERMION-FERMION INTERACTION 

We consider now a more realistic scenario, in which the 
fermionic atoms interacting via a short-range Gaussian 
potential, V(r) = V exp(-&r 2 ) . The mean-field approx- 
imation for the ground-state properties of this system 
has been discussed recently in Ref. [22j , where the BEC- 
BCS crossover was studied by fixing the width of the 
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FIG. 4: (Color online) Energy dependence of the acoustic at- 
tenuation rate as the fermionic component goes through the 
BEC-BCS crossover, for the choice of parameters in Ref. [22l [. 
Along the uig-axis, we plot the position of the energy thresh- 
olds, u) q>T = 2E hta . 
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FIG. 5: (Color online) Effect of the bosonic parameters on 
intensity of the BEC signal. 

potential (r) = 2/v7T& and varying the potential depth 
Vq, or implicitly the scattering length ap of the interac- 
tion. Following Leggett [20| . the dilute fermionic state 
is characterized by the parameter r\ = (kpap)^ 1 , which 
combines together the density of the fermionic gas, char- 
acterized by the Fermi momentum, kp , and the fermion- 
fermion scattering length. The BEC-BCS crossover is 
indicated by the disappearance of the singularity in the 
fermionic density of states, N(k) ~ k 2 | | (see Fig. 3 
in Ref. |22j). These results are similar to a magnetic field 
dependent contact interaction [28|, and are relevant to 
present experiments in which only the lowest hyperfine 
levels of the fermionic atom gas are populated. 

To study the effect of the energy-momentum de- 
pendence of the normal and anomalous densities on 
the acoustic attenuation rate, we consider the BEC- 
BCS crossover at constant-density corresponding to 
kp (r) ps 0.37. The BEC-BCS crossover in this case 
occurs for rj = (kpap)" 1 between 0.15 and 0.80. We 
assume that the dispersion relation of the bosonic com- 
ponent is given by the parameters: c = 0.02 /im MHz, 
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and ttib£,b = 1.6 (u)(/^m). Our results are illustrated 
in Fig. 3J Here, the energy dependence of the atten- 
uation rate, T, as a function of the scattering length 
(or i] = (kpciF) 1 ), is plotted at fixed fermion den- 
sity (kp(r) as 0.37). For our choice of parameters 
in the dispersion relation of the dilute Bose condensate 
(see Eq. @), the ratio q 2 /to q approaches a constant for 
Lj q /eF <C 1. Here kp is the Fermi momentum, and ap is 
the s-wave scattering length in the fermionic component. 

At the threshold, on the BCS side, the acoustic atten- 
uation rate features a sharp edge, followed by a (linear) 
increase in the acoustic attenuation rate as a function 
of the phonon energy; in contrast, on the BEC side, the 
sharp edge is smoothed out, and the acoustic attenua- 
tion rate decreases slowly as a function of energy. It is 
important to recall that in the dilute limit the BEC-BCS 
crossover coincides with the singularity in the scattering 
length, and occurs for fx = [20|. The acoustic atten- 
uation rate provides a crisp signature of the BEC-BCS 
crossover. 

On the BEC side, the energy threshold in first-order 
perturbation theory is higher than the energy threshold 
(u)q t T — 2-Bfe m ) given by energy conservation, due to the 
reduced phase space availability for momentum conserva- 
tion. In this regime we have k m = and the contributions 
due to the third and fourth branches (k < k m ) vanish. 
Due to the lack of available phase space for the inte- 
gral |[7|), the actual threshold in first-order perturbation 
theory is located further away from the energy threshold 
as we go deeper into the BEC regime (see Fig. 0]). The 
effect is enhanced by a smaller bosonic coherence length. 

From the numerical results depicted in Fig. [4] it is ap- 
parent that the attenuation rate calculated in first-order 
perturbation theory vanishes in the deep BEC region. In 
this regime, one needs to include a second-order approxi- 
mation in order to describe the scattering of the boson off 
the molecular condensate. Provided that i) the molecu- 
lar condensate degrees of freedom do not play a role and 
that ii) the coupling between the bosonic and fermionic 
components is sufficiently weak this second-order correc- 
tion in perturbation theory is expected to be small and 
will not be discussed here. For energies greater than the 
threshold, the acoustic attenuation rate becomes nonzero 
and the first-order result is expected to be dominant. 

The dependence on the phonon velocity and/or the 
boson mass of the strength of the BEC signal is illus- 
trated in Fig. [5] Here, we depict the acoustic attenu- 
ation rate for 77 = 0.59 (the crossover region), as ob- 
tained by first changing the phonon velocity from 0.02 to 
0.1 /imMHz, which is equivalent to changing to_b£b from 
1.6 to 0.32 (u)(/Ltm). Next, we increase the boson mass, 
rriB, first by 50 % and then by 100 %, which is equiva- 



lent to reducing by the same factor the bosonic coherence 
length, £b- The acoustic attenuation rate changes little 
when we increase the phonon velocity, c, for a given bo- 
son mass, tub- In contrast, for heavier bosons, the signal 
drops in intensity very quickly for a relatively small in- 
crease of the boson mass, or else for a relatively small 
decrease in the bosonic coherence length. 

We note that in a related formalism, similar shapes (es- 
pecially on the BEC side) are present in the discussion of 
the spin-spin correlation functions (see Refs. [13, H^|). 
These calculations refer to observables accessible via 
Bragg, rf and spin-noise spectroscopy. Therefore, mea- 
surements of the acoustic attenuation rate can be used 
in conjunction with these other techniques to probe fea- 
tures of the fermionic normal and anomalous densities. 
The acoustic attenuation rate discussed here is intrinsi- 
cally a nonzero momentum transfer measurement, thus 
probing different moments of the fermionic normal and 
anomalous densities. This in turn results in different fea- 
tures in the shapes on the BCS side. 

V. CONCLUSIONS 

To summarize, in this paper we have studied the acous- 
tic attenuation rate in the framework of the Fermi-Bosc 
model, for a realistic finite-range fermion-fermion inter- 
action. Our findings suggest that by inducing a density 
fluctuation in the Bose gas component one can induce an 
acoustic response in the fermionic component, and the 
shape of the acoustic attenuation rate can be used as a 
signature of superfluidity on the BCS side. For a low 
excitation energy, the acoustic attenuation rate increases 
with the phonon energy on the BCS side and decreases on 
the BEC side. Therefore, a linear increase of the acous- 
tic attenuation rate with the phonon energy above the 
threshold indicates that one has passed to the BCS side. 
Moreover, our results show that the signal on the BCS 
side is much larger than on the BEC side, which makes 
the detection easier. The acoustic attenuation rate mea- 
surement is favored by a larger phonon velocity and the 
choice of the lightest available boson specie. 
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